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■ It is well known that n integers in the range [Ij^*^] can be sorted in 0{n) time in the 
^ I RAM model using radix sorting. More generally, integers in any range [1, U] can be sorted in 

0{n^/log \ogn) time 5 . However, these algorithms use 0{n) words of extra memory. Is this 
' necessary? 

We present a simple, stable, integer sorting algorithm for words of size O(logn), which works 
in 0{n) time and uses only 0(1) words of extra memory on a RAM model. This is the integer 

■ sorting case most useful in practice. We extend this result with same bounds to the case when 
I the keys are read-only, which is of theoretical interest. Another interesting question is the case 

^ ■ of arbitrary c. Here we present a black-box transformation from any RAM sorting algorithm 

^1 to a sorting algorithm which uses only 0(1) extra space and has the same running time. This 

settles the complexity of in-place sorting in terms of the complexity of sorting. 

■ 1 Introduction 

O 

Given n integer keys S[l . . .n] each in the range [1, n], they can be sorted in 0(n) time using 0{n) 
space by bucket sorting. This can be extended to the case when the keys are in the range 
for some positive constant c by radix sorting that uses repeated bucket sorting with 0{n) ranged 
^ I keys. The crucial point is to do each bucket sorting stably, that is, if positions i and j, i < j, had 

' the same key k, then the copy of k from position i appears before that from position j in the final 

sorted order. Radix sorting takes O (en) time and 0(n) space. More generally, RAM sorting with 
^ \ integers in the range [1, U] is a much-studied problem. Currently, the best known bound is the 

^ • randomized algorithm in [5j that takes 0(n-v/IogTogn) time, and the deterministic algorithm in [T] 

^ that takes O(nloglogn) time. These algorithms also use 0(n) words of extra memory in addition 

to the input. 

We ask a basic question: do we need 0{n) auxiliary space for integer sorting! The ultimate goal 
would be to design in-place algorithms for integer sorting that uses only 0(1) extra words. This 
question has been explored in depth for comparison-based sorting, and after a series of papers, we 
now know that in-place, stable comparison-based sorting can be done in 0(n log n) time [10]. Some 
very nice algorithmic techniques have been developed in this quest. However, no such results are 
known for the integer sorting case. Integer sorting is used as a subroutine in a number of algorithms 
that deal with trees and graphs, including, in particular, sorting the transitions of a finite state 
machine. Indeed, the problem arose in that context for us. In these applications, it is useful if one 
can sort in-place in 0(n) time. From a theoretical perspective, it is likewise interesting to know if 
the progress in RAM sorting, including [31 [U [5], really needs extra space. 

Our results are in-place algorithms for integer sorting. Taken together, these results solve much 
of the issues with space efficiency of integer sorting problems. In particular, our contributions are 
threefold. 
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A practical algorithm. In Section [2l we present a stable integer sorting algorithm for O(logn) 
sized words that takes 0(n) time and uses only 0(1) extra words. 

This algorithm is a simple and practical replacement to radix sort. In the numerous applications 
where radix sorting is used, this algorithm can be used to improve the space usage from 0{n) to 
only 0(1) extra words. We have implemented the algorithm with positive results. 

One key idea of the algorithm is to compress a portion of the input, modifying the keys. The 
space thus made free is used as extra space for sorting the remainder of the input. 

Read-only keys. It is theoretically interesting if integer sorting can be performed in-place without 
modifying the keys. The algorithm above does not satisfy this constraint. In Section [3l we present 
a more sophisticated algorithm that still takes linear time and uses only 0(1) extra words without 
modifying the keys. In contrast to the previous algorithm, we cannot create space for ourselves 
by compressing keys. Instead, we introduce a new technique of pseudo pointers which we believe 
will find applications in other succinct data structure problems. The technique is based on keeping 
a set of distinct keys as a pool of preset read-only pointers in order to maintain linked lists as in 
bucket sorting. 

As a theoretical exercise, in Section \5\ we also consider the case when this sorting has to be 
done stably. We present an algorithm with identical performance that is also stable. Similar to the 
other in-place stable sorting algorithms e.g., comparison-based sorting [10], this algorithm is quite 
detailed and needs very careful management of keys as they are permuted. The resulting algorithm 
is likely not of practical value, but it is still fundamentally important to know that bucket and 
radix sorting can indeed be solved stably in 0{n) time with only 0(1) words of extra space. For 
example, even though comparison-based sorting has been well studied at least since 60's, it was not 
until much later that optimal, stable in-place comparison-based sorting was developed |10j . 

Arbitrary word length. Another question of fundamental theoretical interest is whether the 
recently discovered integer sorting algorithms that work with long keys and sort in o{n log n) time, 
such as [31 [H [5], need any auxiliary space. In Section [U we present a black-box transformation 
from any RAM sorting algorithm to an sorting algorithm which uses only 0(1) extra space, and 
retains the same time bounds. As a result, the running time bounds of [H [5] can now be matched 
with only 0(1) extra space. This transformation relies on a fairly natural technique of compressing 
a portion of the input to make space for simulating space-inefficient RAM sorting algorithms. 

Definitions. Formally, we are given a sequence 5" of n elements. The problem is to sort S 
according to the integer keys, under the following assumptions: 

(i) Each element has an integer key within the interval [1, f/]. 

(ii) The following unit-cost operations are allowed on S: (a) indirect address of any position of 
S; (b) read-only access to the key of any element; (c) exchange of the positions of any two 
elements. 

(iii) The following unit-cost operations are allowed on integer values of 0{\ogU) bits: addition, 
subtraction, bitwise AND/OR and unrestricted bit shift. 

(iv) Only 0(1) auxiliary words of memory are allowed; each word had logC/ bits. 

For the sake of presentation, we will refer to the elements' keys as if they were the input elements. 
For example, for any two elements x, y, instead of writing that the key of x is less than the key 
of y we will simply write x < y. We also need a precise definition of the rank of an element in a 
sequence when multiple occurrences of keys are allowed: the rank of an element Xi in a sequence 
xi . . . xt is the cardinality of the multiset {xj \ xj < xi or [xj = Xi and j < i)}. 
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2 Stable Sorting for Modifiable Keys 

We now describe our simple algorithm for (stable) radix sort without additional memory. 

Gaining space. The first observation is that numbers in sorted order have less entropy than in 
arbitrary order. In particular, n numbers from a universe of u have binary entropy nlogu when 
the order is unspecified, but only log (^) = nlogn — 0(nlogn) in sorted order. This suggests that 
we can "compress" sorted numbers to gain more space: 

Lemma 1. A list of n integers in sorted order can be represented as: (a) an array A[l . . .n] with 
the integers in order; (h) an array of n integers, such that the last (nlogn) hits of the array are 
zero. Furthermore, there exist in-place 0(n) time algorithms for switching between representations 
(a) and (h). 

Proof. One can imagine many representations (b) for which the lemma is true. We note nonetheless 
that some care is needed, as some obvious representations will in fact not lead to in-place encoding. 
Take for instance the appealing approach of replacing A[i] by A[i] — S[i — 1], which makes numbers 
tend to be small (the average value is ^). Then, one can try to encode the difference using a code 
optimized for smaller integers, for example one that represents a value x using logx + O(loglogx) 
bits. However, the obvious encoding algorithm will not be in-place: even though the scheme is 
guaranteed to save space over the entire array, it is possible for many large values to cluster at the 
beginning, leading to a rather large prefix being in fact expanded. This makes it hard to construct 
the encoding in the same space as the original numbers, since we need to shift a lot of data to the 
right before we start seeing a space saving. 

As it will turn out, the practical performance of our radix sort is rather insensitive to the exact 
space saving achieved here. Thus, we aim for a representation which makes in-place encoding 
particularly easy to implement, sacrificing constant factors in the space saving. 

First consider the most significant bit of all integers. Observe that if we only remember the 
minimum i such that A[i] > m/2, we know all most significant bits (they are zero up to i and one 
after that). We will encode the last n/3 values in the array more compactly, and use the most 
significant bits of ^[1 . . . |n] to store a stream of |n bits needed by the encoding. 

We now break a number x into hi(j;), containing the upper [log2(?^/3)J bits, and lo(x), with 
the low logn — [log2(n/3)J bits. For all values in ^[|n -|- 1 . . .n], we can throw away hi(j4[z]) as 
follows. First we add hi(A[|n -|- 1]) zeros to the bit stream, followed by a one; then for every 
i = |n -|- 2, . . . ,n we add hi(A[z]) — hi(^[i — 1]) zeros, followed by a one. In total, the stream 
contains exactly n/3 ones (one per element), and exactly hi(^[n]) < n/3 zeros. Now we simply 
compact lo(A[|n -|- 1]), . . . , lo(A[n]) in one pass, gaining | [log2(n/3)J free bits. □ 

An unstable algorithm. Even just this compression observation is enough to give a simple 
algorithm, whose only disadvantage is that it is unstable. The algorithm has the following structure: 

1. sort the subsequence . . . (n/logn)] using the optimal in-place mergesort in jlOj . 

2. compress S\i . . . (n/logn)] by Lemma[Tl generating Vt{n) bits of free space. 

3. radix sort S[{n/ log n) -|- 1 . . . n] using the free space. 

4. uncompress 5 [1 ... (n/ log n)]. 

5. merge the two sorted sequences S\i . . . (nj logn)] and 5'[(n/ logn) -|- 1 . . . n] by using the in- 
place, linear time merge in |10j . 



3 



The only problematic step is 3. The implementation of this step is based on the cycle leader 
permuting approach where a sequence A is re-arranged by following the cycles of a permutation 
TT. First A[l] is sent in its final position vr(l). Then, the element that was in vr(l) is sent to its 
final position 7r(7r(l)). The process proceeds in this way until the cycle is closed, that is until the 
element that is moved in position 1 is found. At this point, the elements starting from A[2] are 
scanned until a new cycle leader A[i] (i.e. its cycle has not been walked through) is found, ^[i]'s 
cycle is followed in its turn, and so forth. 

To sort, we use 2n'^ counters ci, . . . , c^e and di, . . . , d„e . They are stored in the auxiliary words 
obtained in step 2. Each dj is initialized to 0. With a first scan of the elements, we store in any q 
the number of occurrences of key i. Then, for each i = 2 . . . rf, we set Cj = Cj_i + 1 and finally we 
set ci = 1 (in the end, for any i we have that Cj = X]j<j + Now we have all the information 
for the cycle leader process. Letting j = (n/logn) + 1, we proceed as follows: 

(i) let i be the key of S[j]; 

(a) if Ci < j < Ci+i then S[j] is already in its final position, hence we increment j by 1 and go to 
step (i); 

{iii) otherwise, we exchange S[j] with S[ci + di], we increment di by 1 and we go to step (i). 

Note that this algorithm is inherently unstable, because we cannot differentiate elements which 
should fall between Cj and q+i — 1, given the free space we have. 

Stability through recursion. To achieve stability, we need more than n free bits, which we can 
achieve by bootstrapping with our own sorting algorithm, instead of merge sort. There is also an 
important practical advantage to the new stable approach: the elements are permuted much more 
conservatively, resulting in better cache performance. 

1. recursively sort a constant fraction of the array, say S[l . . . n/2]. 

2. compress ^[l . . . n/2] by LemmalU generating r2(nlogn) bits of free space. 

3. for a small enough constant 7, break the remaining n/2 elements into chunks of jn numbers. 
Each chunk is sorted by a classic radix sort algorithm which uses the available space. 

4. uncompress ^[l ... n/2]. 

5. we now have 1 + I/7 = 0(1) sorted subarrays. We merge them in linear time using the stable 
in-place algorithm of [lOj. 

We note that the recursion can in fact be implemented bottom up, so there is no need for a stack 
of superconstant space. For the base case, we can use bubble sort when we are down to n < ^/no 
elements, where no is the original size of the array at the top level of the recursion. 

Steps 2 and 4 are known to take 0{n) time. For step 3, note that radix sort in base R applied 
to numbers requires N + R additional words of space, and takes time 0{Nlogj^u). Since we 
have a free space of r2(nlogn) bits or 0,(n) words, we can set N = R = ^n, for a small enough 
constant 7. As we always have n = ^(y'no) = u^^^\ radix sort will take linear time. 

The running time is described by the recursion T(n) = T(n/2) + 0(n), yielding T(n) = 0(n). 
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A self-contained algorithm. Unfortunately, all algorithms so far use in-place stable merging 
algorithm as in [10] . We want to remove this dependence, and obtain a simple and practical sorting 
algorithm. By creating free space through compression at the right times, we can instead use a 
simple merging implementation that needs additional space. We first observe the following: 

Lemma 2. Let k >2 and a > be arbitrary constants. Given k sorted lists of n/k elements, and 
an words of free space, we can merge the lists in 0{n) time. 

Proof. We divide space into blocks of an/{k + 1) words. Initially, we have A; + 1 free blocks. We 
start merging the lists, writing the output in these blocks. Whenever we are out of free blocks, 
we look for additional blocks which have become free in the original sorted lists. In each list, the 
merging pointer may be inside some block, making it yet unavailable. However, we can only have 
k such partially consumed blocks, accounting for less than k^^ wasted words of space. Since in 
total there are an free words, there must always be at least one block which is available, and we 
can direct further output into it. 

At the end, we have the merging of the lists, but the output appears in a nontrivial order of 
the blocks. Since there are {k + 1)(1 + 1/a) = 0(1) blocks in total, we can remember this order 
using constant additional space. Then, we can permute the blocks in linear time, obtaining the 
true sorted order. □ 

Since we need additional space for merging, we can never work with the entire array at the 
same time. However, we can now use a classic sorting idea, which is often used in introductory 
algorithms courses to illustrate recursion (see, e.g. [2]). To sort n numbers, one can first sort the 
first |n numbers (recursively), then the last |n numbers, and then the first |n numbers again. 
Though normally this algorithm gives a running time of a;(n^), it works efficiently in our case 
because we do not need recursion: 

1. sort ^[l . . . n/3] recursively. 

2. compress ^[l •••§], and sort ^[l^ + l ... n] as before: first radix sort chunks of 771 numbers, and 
then merge all chunks by Lemma [2] using the available space. Finally, uncompress S[l •••§]• 

3. compress + 1 . . . n], which is now sorted. Using Lemma [21 merge ^[l • • • §] with + 
1 . . . Finally uncompress. 

4. once again, compress S[l •••§], merge + 1 . . . ^] with + 1 . . . n], and uncompress. 

Note that steps 2-4 are linear time. Then, we have the recursion T(n) = T(n/3) + 0{n), solving 
to T{n) = 0{n). Finally, we note that stability of the algorithm follows immediately from stability 
of classic radix sort and stability of merging. 

Practical experience. The algorithm is surprisingly effective in practice. It can be implemented 
in about 150 lines of C code. Experiments with sorting 1-10 million 32-bit numbers on a Pentium 
machine indicate the algorithm is roughly 2.5 times slower than radix sort with additional memory, 
and slightly faster than quicksort (which is not even stable). 

3 Unstable Sorting for Read-only Keys 
3.1 Simulating auxiliary bits 

With the bit stealing technique [9j, a bit of information is encoded in the relative order of a pair 
of elements with different keys: the pair is maintained in increasing order to encode a and vice 
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versa. The obvious drawback of this technique is that the cost of accessing a word of w encoded 
bits is O (w) in the worst case (no word-level parallelism). However, if we modify an encoded word 
with a series of / increments (or decrements) by 1, the total cost of the entire series is O (l) (see pj). 

To find pairs of distinct elements, we go from S* to a sequence Z'Y'XY"Z" with two properties. 
(i) For any z' G Z', y' £ Y' , x G X, y" £ Y" and z" £ Z" we have that z' < y' < x < y" < z" . (ii) 
Let m = a [n/ log n] , for a suitable constant a. Y' is composed by the element y'^ with rank m plus 
all the other elements equal to y'^. Y" is composed by the element y'^^ with rank n—m+1 plus all the 
other elements equal to y'^. To obtain the new sequence we use the in-place, linear time selection 
and partitioning algorithms in [6l[7j. If X is empty, the task left is to sort Z' and Z" , which can be 
accomplished with any optimal, in-place mergesort (e.g. [10]. Let us denote Z'Y' with M' and Y" Z" 
with M". The m pairs of distinct elements (M'[l], M"[l]), (M'[2], M"[2]), . . . , {M'[m], M"[m]) wih 
be used to encode information. 

Since the choice of the constant a does not affect the asymptotic complexity of the algorithm, 
we have reduced our problem to a problem in which we are allowed to use a special hit memory 
with 0(re/logn) bits where each bit can be accessed and modified in constant time but without 
word-level parallelism. 



3.2 Simulating auxiliary memory for permuting 

With the internal buffering technique [8j, some of the elements are used as placeholders in order 
to simulate a working area and permute the other elements at lower cost. In our unstable sorting 
algorithm we use the basic idea of internal buffering in the following way. Using the selection and 
partitioning algorithms in [6l[7j, we pass from the original sequence S to ABC with two properties. 
[i) For any a £ h £ B and c G C, we have that a <h < c. {ii) B is composed of the element b' 
with rank \n/2\ plus all the other elements equal to b' . We can use BC as an auxiliary memory 
in the following way. The element in the first position of BC is the separator element and will 
not be moved. The elements in the other positions of BC are placeholders and will be exchanged 
with (instead of being overwritten by) elements from A in any way the computation on A (in our 
case the sorting of A) may require. The "emptiness" of any location i of the simulated working 
area in BC can be tested in 0(1) time by comparing the separator element -BC[1] with BC[i]: if 
-BC[1] < BC[i] the ith location is "empty" (that is, it contains a placeholder), otherwise it contains 
one of the elements in A. 

Let us suppose we can sort the elements in A in 0(1^41) time using BC as working area. After 
A is sorted we use the partitioning algorithm in [6] to separate the elements equal to the separator 
element (-BC[1]) from the elements greater than it (the computation on A may have altered the 
original order in BC). Then we just re-apply the same process to C, that is we divide it into 
A'B'C, we sort A' using B'C as working area and so forth. Clearly, this process requires 0{n) 
time and when it terminates the elements are sorted. Obviously, we can divide A into p = O {1) 
equally sized subsequences Ai, A2 ■ ■ ■ Ap, then sort each one of them using BC as working area and 
finally fuse them using the in-place, linear time merging algorithm in [10]. Since the choice of the 
constant p does not affect the asymptotic complexity of the whole process, we have reduced our 
problem to a new problem, in which we are allowed to use a special exchange memory of O (n) 
locations, where each location can contain input elements only (no integers or any other kind of 
data). Any element can be moved to and from any location of the exchange memory in 0(1) time. 



3.3 The reduced problem 

By blending together the basic techniques seen above, we can focus on a reduced problem in which 
assumption (iv) is replaced by: 
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(iv) Only 0(1) words of normal auxiliary memory and two kinds of special auxiliary memory are 
allowed: 



(a) A random access bit memory ^ with O (n/logn) bits, where each bit can be accessed in 
0(1) time (no word-level parallelism). 

(b) A random access exchange memory S with O {n) locations, where each location can 
contain only elements from S and they can be moved to and from any location of S in 
0(1) time. 

If we can solve the reduced problem in O (n) time we can also solve the original problem with the 
same asymptotic complexity. However, the resulting algorithm will be unstable because of the use 
of the internal buffering technique with a large pool of placeholder elements. 

3.4 The naive approach 

Despite the two special auxiliary memories, solving the reduced problem is not easy. Let us consider 
the following naive approach. We proceed as in the normal bucket sorting: one bucket for each 
one of the range values. Each bucket is a linked list: the input elements of each bucket are 
maintained in S while its auxiliary data (e.g. the pointers of the list) are maintained in SS. In 
order to amortize the cost of updating the auxiliary data (each pointer requires a word of (log n) 
bits and does not have word- level parallelism) , each bucket is a linked list of slabs of (log^ n) 
elements each (,^ has only O (n/logn) bits). At any time each bucket has a partially full head slab 
which is where any new element of the bucket is stored. Hence, for each bucket we need to store 
in =^ a word of O (log log n) bits with the position in the head slab of the last element added. The 
algorithm proceeds as usual: each element in S is sent to its bucket in O (1) time and is inserted in 
the bucket's head slab. With no word-level parallelism in ^ the insertion in the head slab requires 
O (log log n) time. Therefore, we have an 0(n log log n) time solution for the reduced problem and, 
consequently, an unstable 0(n log log n) time solution for the original problem. 

This simple strategy can be improved by dividing the head slab of a bucket into second level 
slabs of Q (log log n) elements each. As for the first level slabs, there is a partially full, second level 
head slab. For any bucket we maintain two words in ^: the first one has O (log log log n) bits and 
stores the position of the last element inserted in the second level head slab; the second one has 
O (log log n) bits and stores the position of the last full slab of second level contained in the first 
level head slab. Clearly, this gives us an O(nlogloglogn) time solution for the reduced problem 
and the corresponding unstable solution for the original problem. By generalizing this approach to 
the extreme, we end up with O (log* n) levels of slabs, an O (n log* n) time solution for the reduced 
problem and the related unstable solution for the original problem. 

3.5 The pseudo pointers 

Unlike bit stealing and internal buffering which were known earlier, the pseudo pointers technique 
has been specifically designed for improving the space complexity in integer sorting problems. 
Basically, in this technique a set of elements with distinct keys is used as a pool of pre-set, read- 
only pointers in order to simulate efficiently traversable and updatable linked lists. Let us show 
how to use this basic idea in a particular procedure that will be at the core of our optimal solution 
for the reduced problem. 

Let d be the number of distinct keys in S. We are given two sets of d input elements with 
distinct keys: the sets Q and V of guides and pseudo pointers, respectively. The guides are given 
us in sorted order while the pseudo pointers form a sequence in arbitrary order. Finally, we are 
given a multiset X of d input elements (i.e. two elements of I can have equal keys). The procedure 
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uses the guides, the pseudo pointers and the exchange memory to sort the d input elements of X in 
O (d) time. 

We use three groups of contiguous locations in the exchange memory The first group H 
has n*^ locations (one for each possible value of the keys). The second group L has n"^ slots of two 
adjacent locations each. The last group R has d locations, the elements of I will end up here in 
sorted order. H, L and R are initially empty. We have two main steps. 

First. For each s G T, we proceed as follows. Let p be the leftmost pseudo pointer still in V. 
If the sth location of H is empty, we move p from V to H[s] and then we move s from I to the 
first location of L\p\ (i.e. the first location of the pth slot of L) leaving the second location of L\p] 
empty. Otherwise, if H[s] contains an element p' (a pseudo pointer) we move s from I to the first 
location of L\p], then we move p' from H[s] to the second location of L[p] and finally we move p 
from V to H[s]. 

Second. We scan the guides in Q from the smallest to the largest one. For a guide g G we 
proceed as follows. If the gth location of H is empty then there does not exist any element equal to 
g among the ones to be sorted (and initially in 2) and hence we move to the next guide. Otherwise, 
if H[G] contains a pseudo pointer p, there is at least one element equal to g among the ones to 
be sorted and this element is currently stored in the first location of the pth slot of L. Hence, we 
move that element from the first location of L[p] to the leftmost empty location of R. After that, 
if the second location of L\p] contains a pseudo pointer p', there is another element equal to g and 
we proceed in the same fashion. Otherwise, if the second location of L[p] is empty then there are 
no more elements equal to g among the ones to be sorted and therefore we can focus on the next 
guide element. 

Basically, the procedure is bucket sorting where the auxiliary data of the list associated to each 
bucket (i.e. the links among elements in the list) is implemented by pseudo pointers in V instead 
of storing it explicitly in the bit memory (which lacks of word-level parallelism and is inefficient 
in access). It is worth noting that the buckets' lists implemented with pseudo pointers are spread 
over an area that is larger than the one we would obtain with explicit pointers (that is because 
each pseudo pointer has a key of logn'^ bits while an explicit pointer would have only logd bits). 

3.6 The optimal solution 

We can now describe the algorithm, which has three main steps. 

First. Let us assume that for any element s G S there is at least another element with the same 
key. (Otherwise, we can easily reduce to this case in linear time: we isolate the O [n'') elements 
that do not respect the property, we sort them with the in-place mergesort in [lOj and finally we 
merge them after the other 0{n) elements are sorted.) With this assumption, we extract from S 
two sets Q and V of d input elements with distinct keys (this can be easily achieved in O (n) time 
using only the exchange memory S") . Finally we sort Q with the optimal in-place mergesort in [TO] . 

Second. Let S' be the sequence with the (O (n)) input elements left after the first step. Using 
the procedure in § 13.51 (clearly, the elements in the sets Q and V computed in the first step will be 
the guides and pseudo pointers used in the procedure), we sort each block Bi of S' with d contiguous 
elements. After that, let us focus on the first t = O (log log n) consecutive blocks Bi, B2, . . . , Bt. 
We distribute the elements of these blocks into < t groups Gi , G2 . . . in the following way. Each 
group Gj can contain between d and 2d elements and is allocated in the exchange memory The 
largest element in a group is its pivot. The number of elements in a group is stored in a word of 
Q (log d) bits allocated in the bit memory Initially there is only one group and is empty. In the 
ith. step of the distribution we scan the elements of the ith. block Bi . As long as the elements of Bi 
are less than or equal to the pivot of the first group we move them into it. If, during the process, 
the group becomes full, we select its median element and partition the group into two new groups 
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(using the selection and partitioning algorithms in [6l[7I). When, during the scan, the elements of 
Bi become greater than the pivot of the first group, we move to the second group and continue in 
the same fashion. It is important to notice that the number of elements in a group (stored in a 
word of (log d) bits in the bit memory ^) is updated by increments by 1 (and hence the total 
cost of updating the number of elements in any group is linear in the final number of elements 
in that group, see [2]). Finally, when all the elements of the first t = (log log n) consecutive 
blocks Bi, B2, ■ ■ ■ , Bt have been distributed into groups, we sort each group using the procedure 
in § 13.51 (when a group has more than d elements, we sort them in two batches and then merge 
them with the in-place, linear time merging in |10j). The whole process is repeated for the second 
t = Q (log log n) consecutive blocks, and so forth. 

Third. After the second step, the sequence S' (which contains all the elements of S with 
the exclusion of the guides and pseudo pointers, see the first step) is composed by contiguous 
subsequences S[,S2, - ■ ■ which are sorted and contain G(dloglogn) elements each (where d is the 
number of distinct elements in S). Hence, if we see S' as composed by contiguous runs of elements 
with the same key, we can conclude that the number of runs of S' is O (n/ log log n). Therefore 
S' can be sorted in O (n) time using the naive approach described in § 13.41 with only the following 
simple modification. As long as we are inserting the elements of a single run in a bucket, we 
maintain the position of the last element inserted in the head slab of the bucket in a word of 
auxiliary memory (we can use 0(1) of them) instead of accessing the inefficient bit memory ^ at 
any single insertion. When the current run is finally exhausted, we copy the position in the bit 
memory. Finally, we sort V and we merge V, A and S' (once again, using the sorting and merging 
algorithms in [TO]). 

3.7 Discussion: Stability and Read-only Keys 

Let us focus on the reasons why the algorithm of this section is not stable. The major cause of insta- 
bility is the use of the basic internal buffering technique in conjunction with large {uj{ polylog (n))) 
pools of placeholder elements. This is clearly visible even in the first iteration of the process in 
§ 13.21 after being used to permute A into sorted order, the placeholder elements in BC are left 
permuted in a completely arbitrary way and their initial order is lost. 

4 Reducing Space in any RAM Sorting Algorithm 

In this section, we consider the case of sorting integers of w = a;(logn) bits. We show a black 
box transformation from any sorting algorithm on the RAM to a stable sorting algorithm with the 
same time bounds which only uses 0(1) words of additional space. Our reduction needs to modify 
keys. Furthermore, it requires randomization for large values of w. 

We first remark that an algorithm that runs in time t{n) can only use 0{t{n)) words of space 
in most realistic models of computation. In models where the algorithm is allowed to write t{n) 
arbitrary words in a larger memory space, the space can also be reduced to 0{t{n)) by introducing 
randomization, and storing the memory cells in a hash table. 

Small word size. We first deal with the case w = polylog (n). The algorithm has the following 
structure: 

1. sort S[\ . . . n/ log n] using in-place stable merge sort [lOj . Compress these elements by Lemma 
[1] gaining Q{n) bits of space. 
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2. since = 0(n log n), the RAM sorting algorithm uses at most 0(t(n)-ri;) = 0(n polylog (n)) 
bits of space. Then we can break the array into chunks of n/log'^n elements, and sort each 
one using the available space. 

3. merge the log'^ n sorted subarrays. 

4. uncompress S[l . . . n/ logn] and merge with the rest of the array by stable in-place merg- 
ing [TO]. 

Steps 1 and 4 take linear time. Step 2 requires log'^ n • t(n/ log"^ n) = 0{t{n)) because t{n) is convex 
and bounded in [n,nlogn]. We note that step 2 can always be made stable, since we can afford a 
label of O(logn) bits per value. 

It remains to show that step 3 can be implemented in 0{n) time. In fact, this is a combination 
of the merging technique from Lemma [2] with an atomic heap [1] . The atomic heap can maintain 
a priority queue over polylog (n) elements with constant time per insert and extract-min. Thus, 
we can merge log"^ n lists with constant time per element. The atomic heap can be made stable by 
adding a label of clog log n bits for each element in the heap, which we have space for. The merging 
of Lemma [2] requires that we keep track of 0{k/a) subarrays, where k = log'^n was the number of 
lists and a = 1/ polylog (n) is fraction of additional space we have available. Fortunately, this is 
only polylog (n) values to record, which we can afford. 

Large word size. For word size w > log^'^^^n, the randomized algorithm of [Ij can sort in 0{n) 
time. Since this is the best bound one can hope for, it suffices to make this particular algorithm 
in-place, rather than give a black-box transformation. We use the same algorithm from above. The 
only challenge is to make step 2 work: sort n keys with 0(n polylog (n)) space, even if the keys 
have w > polylog (n) bits. 

We may assume w > log^ n, which simplifies the algorithm of [Ij to two stages. In the first 
stage, a signature of O(log^n) bits is generated for each input value (through hashing), and these 
signatures are sorted in linear time. Since we are working with 0(log'^ n)-bit keys regardless of the 
original w, this part needs 0(n polylog (n)) bits of space, and it can be handled as above. 

From the sorted signatures, an additional pass extracts a subkey of t(;/logn bits from each 
input value. Then, these subkeys are sorted in linear time. Finally, the order of the original keys 
is determined from the sorted subkeys and the sorted signatures. 

To reduce the space in this stage, we first note that the algorithm for extracting subkeys does 
not require additional space. We can then isolate the subkey from the rest of the key, using shifts, 
and group subkeys and the remainder of each key in separate arrays, taking linear time. This way, 
by extracting the subkeys instead of copying them we require no extra space. We now note that 
the algorithm in [1] for sorting the subkeys also does not require additional space. At the end, we 
recompose the keys by applying the inverse permutation to the subkeys, and shifting them back 
into the keys. 

Finally, sorting the original keys only requires knowledge of the signatures and order information 
about the subkeys. Thus, it requires 0(n polylog (n)) bits of space, which we have. At the end, we 
find the sorted order of the original keys and we can implement the permutation in linear time. 

5 Stable sorting for read-only keys 

In the following we will denote n'' with r. Before we begin, let us recall that two consecutive 
sequences X and Y, possibly of different sizes, can be exchanged stably, in-place and in linear 
time with three sequence reversals, since YX = (X^Y^)^. Let us give a short overview of the 
algorithm. We have three phases. 
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Preliminary phase (§ 15. ip . The purpose of this phase is to obtain some cohections of elements 
to be used with the three techniques described in §[3l We extract (n/logn) smallest and largest 
elements of S. They will form an encoded memory of Q (n/logn) bits. Then, we extract from the 
remaining sequence (n*^) smallest elements and divide them into 0(1) jump zones of equal length. 
After that, we extract from the remaining sequence some equally sized sets of distinct elements. 
Each set is collected into a contiguous zone. At the end of the phase, we have guide, distribution, 
pointer, head and spare zones. 

Aggregating phase (§ 15. 2|) . After the preliminary phase, we have reduced the problem to 
sorting a smaller sequence S' (still of O (n) size) using various sequences built to be used with 
the basic techniques in § [3j Let d be the number of distinct elements in S' (computed during the 
preliminary phase). The objective of this phase is to sort each subsequence of size Q (dpolylog (n)) 
of the main sequence S' . For any such subsequence S'l, we first find a set of pivots and then sort S'l 
with a distributive approach. The guide zone is sorted and is used to retrieve in sorted order lists of 
equal elements produced by the distribution. The distribution zone provides sets of pivots elements 
that are progressively moved into one of the spare zones. The head zone furnishes placeholder 
elements for the distributive processes. The distributive process depends on the hypothesis that 
each d contiguous elements of S'l are sorted. The algorithm for sorting (d) contiguous elements 
stably, in O (d) time and O (1) space (see § 15.2. ip employs the pseudo pointers technique and the 
guide, jump, pointer and spare zones are crucial in this process. 

Final phase (§ 15. Sh . After the aggregating phase the main sequence S' has all its subsequences 
of (dpolylog (n)) elements in sorted order. With an iterative merging process, we obtain from S' 
two new sequences: a small sequence containing O (dlog^ n) sorted elements; a large sequence still 
containing O (n) elements but with an important property: the length of any subsequence of equal 
elements is multiple of a suitable number (log^ n) . By exploiting its property, the large sequence 
is sorted using the encoded memory and merged with the small one. Finally, we take care of all 
the zones built in the preliminary phase. Since they have sizes either O (n/logn) or O (n*^), they 
can be easily sorted within our target bounds. 

5.1 Preliminary Phase 

The preliminary phase has two main steps described in Sections 15.1.11 and l5.1.2[ 
5.1.1 Encoded memory 

We start by collecting some pairs of distinct elements. We go from to a sequence Z'Y' XY" Z" 
with the same two properties we saw in § 13.11 We use the linear time selection and partitioning 
in [6l [7] which are also stable. Let us maintain the same notations used in § 13.11 We end up 
with m pairs of distinct elements (M'[l], M"[l]), . . . , (M' [m], M" [m]) to encode information by bit 
stealing. We use the encoded memory based on these pairs as if it were actual memory (that is, we 
will allocate arrays, we will index and modify entries of these arrays, etc). However, in order not 
to lose track of the costs, the names of encoded structures will always be written in the following 
way: I, U, etc. 

We allocate two arrays If,^ and len, each one with r = entries of 1 bit each. The entries 
of both arrays are set to 0. Ibg and len will be used each time the procedure in § 15.2.11 of the 
aggregating phase is invoked. 
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5.1.2 Jump, guide, distribution, pointer, head and spare zones 

The second main step of this phase has six sub-steps. 

First. Let us suppose \X\ > n/logn (otherwise we sort it with the mergesort in |10]). Using the 
selection and partitioning in [6l [3, we go from X to JX' such that J is composed by the element j* 
with rank 3r + 1 = 3n^ + 1 (in X) plus all the elements (in X) < j* . Then, we move the rightmost 
element equal to j* in the last position of J (easily done in O (| J|) and stably with a sequence 
exchange) . 

Second. Let us suppose \X'\ > n/logn. With this step and the next one we separate the ele- 
ments which appear more than 7 times. Let us allocate in our encoded memory of m = O (n/logn) 
bits an encoded array I with r (= n^) entries of 4 bits each. All the entries are initially set to 0. 
Then, we start scanning X' from left to right. For any element u ^ X' accessed during the scan, if 
l[u] < 7, we increment I[u] by one. 

Third. We scan X' again. Let u £ X' he the ith element accessed. If I[u] < 7, we decrement 
I[u] by 1 and exchange X'[i] (= u) with J[i]. At the end of the scan we have that J = WJ", where 
W contains the elements of X' occurring less than 7 times in X' . Then, we have to gather the 
elements previously in J and now scattered in X' . We accomplish this with the partitioning in [6], 
using J[| J|] to discern between elements previously in J and the ones belonging to X' (we know 
that J[\J\] is equal to j* and, for any j £ J and any x' £ X', j < J[\J\] < x'). After that we have 
W J" J' X" where the elements in J' are the ones previously in J and exchanged during the scan of 
X' . We exchange W with J' ending up with JWX" . 

Fourth. We know that each element of X" occurs at least 7 times in it. We also know that the 
entries of I encode either or 7. We scan X" from left to right. Let u £ X" be the ith element 
accessed. If I[u] = 7, we decrement I[u] by one and we exchange X"[i] (= u) with J[i]. After the 
scan we have that J = GJ'" , where, for any j, G[j] was the leftmost occurrence of its kind in X" 
(before the scan). Then, we sort G with the mergesort in [10] (|G| = O (r) = O (n^) and e < 1). 
After that, similarly to the third step, we gather the elements previously in J and now scattered in 
X" because of the scan. We end up with the sequence JWGX'" . We repeat the same process (only 
testing for I[u] equal to 6 instead of 7) to gather the leftmost occurrence of each distinct element 
in X'" into a zone D, ending up with the sequence JWGDX"" . 

Fifth. Each element of X"" occurs at least 5 times in it and the entries of I encode either or 
5. We scan X"" , let u £ X"" be the ith element accessed. If l[u] = 5, we decrement l[v\ by 1 and 
exchange (= u) with J[i]. After the scan we have that J = PJ'" , where, for any j, P[j] was 

the leftmost occurrence of its kind in X"" (before the scan). Unlike the fourth step, we do not sort 
P. We repeat the process finding Ti, T2, T3 and H containing the second, third, fourth and fifth 
leftmost occurrence of each distinct element in X"" , respectively. After any of these processes, we 
gather back the elements previously in J scattered in X"" (same technique used in the third and 
fourth steps). We end up with the sequence JWGDPT1T2T3HS' . 

Sixth. Let us divide J into J1J2J3V, where | Ji| = | J2I = | J3I = r and \V\ = \ J\ — 3r. We scan 
G, let u G G be the ith element accessed, we exchange Ti[i], T2[i] and T^li] with Ji[u], and 
Jslu], respectively. 

5.1.3 Summing up 

We will refer to Ji, J2 and J3 as jump zones. Zone G, D, P and H will be referred to as guide, distri- 
bution, pointer and head zones, respectively. Finally, Ti, T2 and T3 will be called spare zones. With 
the preliminary phase we have passed from the initial sequence S to M' JiJ2J3VWGDPTiT2T^HS' M' 
We allocated in the encoded memory two arrays \g and Ign- The encoded memory, \g and Ign, and 
the jump, guide, distribution, pointer, head and spare zones will be used in the next two phases 
to sort the S' . Zones V and W are a byproduct of the phase and will not have an active role 
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in the sorting of S'. The number of distinct elements in sequence S' is less than or equal to the 
sizes of guide, distribution, pointer and head zones. For the rest of the paper we will denote |G| 
(= \D\ = \P\ = \H\) with d. 

Lemma 3. The preliminary phase requires 0{n) time, uses 0(1) auxiliary words and is stable. 
5.2 Aggregating Phase 

Let us divide S' into k subsequences 5*2 ... 5^ with \Sl\ = dlog'^ n, for a suitable constant /3 > 4. 
Let t = log^ n, for a suitable constant 6 < 1. We will assume that d> {2t + 1) log \S'j^\ . We leave the 
particular case where d < (2t + 1) log \ S^\ for the full paper. For a generic 1 < / < A;, let us assume 
that any 5"^ with /' < / has been sorted and that H is next to the left end of S'^. To sort S'l we have 
two main steps described in § 15.2.21 and § 15.2.31 They rely on the algorithm described in § 15.2.11 

5.2.1 Sorting O (d) contiguous elements 

Let us show how to exploit the two arrays 1^^ and len, (in the encoded memory in the preliminary 
phase) and the jump, guide and pointer zones to sort a sequence A, with 1^41 < d, stably in O {\A\) 
time and using O (1) auxiliary words. The process has two steps. 

First. We scan A, let u £ A be the ith element accessed. Let p = P[i] and h = Ji[u]. If 
Ibg[u] = 0, we set both Ibg[u] and Ien\p] to 1. In any case, we exchange Ji[u\ (= h) with J2[p] and 
A[i] (= u) with JsIp]. Then, we exchange P[i] {= p) with Ji[u] (which is not h anymore). 

Second. Let j = \A\. We scan G, let 5 G G be the ith element accessed. If Ibg[g] = 0, we do 
nothing. Otherwise, let p be Ji[g], we set Ibgid] = and execute the following three steps, (i) We 
exchange Jslp] with A\j], then Ji[g] with P[j] and finally Ji[g] with J2\p]. [ii) We decrease j by 1. 
{Hi) If len[p] = 1, we set \n\p\ = and the sub-process ends, otherwise, let p be Ji[g + 1], and we 
go to (i). 

Let us remark that the O (|A|) entries of \g and len that are changed from to 1 in the first 
step, are set back to in the second one. 1;,^ and len have been initialized in the preliminary phase. 
We could not afford to re-initialize them every time we invoke the process (they have r = n*^ entries 
and |j4| may be o(n^)). 

Lemma 4. Using the encoded arrays \g, len oind the jump, guide and pointer zones, a sequence A 
with \A\ < d can he sorted stably, in 0{\A\) time and using 0(1) auxiliary words. 

5.2.2 Finding pivots for S'l 

We find a set of pivots {ei, 62 ... , ep_i, ep} with the following properties: (i) |{x S S"; | x < ei}| < 
d; {ii) \{x £ S'i\ Cp < x}\ < d; {Hi) \{x £ S'i\ Ci < x < ei+i}| < d, for any 1 < i < p; {iv) p = 
Q (log'^n). In the end the pivots reside in the first p positions of D. We have four steps. 

First. We allocate in the encoded memory an array P with r entries of log \ S'i\ bits, but we do 
not initialize each entry of P. We initialize the only d of them we will need: for any i = 1 . . . |G|, 
we set ¥[G[i\] to 0. 

Second. We scan S'^ from left to right. Let u G S'^ be the ith element accessed, we increment 
P[n] by 1. 

Third. We sort the distribution zone D using the algorithm described in § 15.2. 1[ 

Fourth. Let i = 1, j = and p = 0. We repeat the following process until i > \G\. (i) Let 

u = G[i], we set j = j + F[u]. {ii) If j < d we increase i by 1 and go to (i). If j > d, we increase p 

by 1, exchange D[i] with D\p], increase i by 1, set j to and go to (i). 
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5.2.3 Sorting S'l 



Let p be the number of pivots for S'^ selected in § 15.2.21 and now residing in the first p positions of 
D. Let u be log Let us assume that H is next to the left end of 5-. We have six steps. 

First. Let i = and j = 0. The following two steps are repeated until i > p: (i) we increase i 
by p/t and j by 1; (ii) we exchange D[i] with Ti[j]. We will denote with p' the number of selected 
pivots, now temporarily residing in the first p' positions of Ti. 

Second. Let us divide S'l into q = \S'i \ /d blocks B1B2 . . . Bq of d elements each. We sort each 
Bi using the algorithm in § I5.2.1[ 

Third. With a sequence exchange we bring H next to the right end of S[. Let us divide H 
into HiHi . . . Hp'HpiHpi^iH' , where = \Hi\ = \Hi\ = u. Let / = /u. We allocate the 

following arrays: (i) U^mc and Upre both with / + 2p' + 1 entries of Q (u) bits; (ii) EI and M with 
p' + 1 and p' entries of Q (logu) bits; {Hi) L and L with p' + 1 and p' entries of Q (u) bits; {iv) N 
and N with p' + 1 and p' entries of G (u) bits. Each entry of any array is initialized to 0. 

Fourth. In this step we want to transform S'l and H in the following ways. We pass from S'l to 
U1U2 ■ ■ ■ Uf'-iUf'H" , where the C/j's are called units, for which the following holds. 

(i) f >f - {2p' + 1) and \Ui\ = u, for any 1 < i < /'. 

(ii) For any Ui, 1 < i < f , one of the following holds: (a) there exists a 1 < j < p' such that 
X = Ti[j], for any x £ Ui; (6) x < Ti[l], for any x £ Ui; (c) Ti\p'] < x, for any x S [/«; (d) 
there exists a 1 < j' < p' — 1 such that Ti[j'] < x < Ti[j' + 1], for any x £ Ui. 

(Hi) Let us call a set of related units a maximal set of units U = {Ui-^ , Ui.^ , • • • , Ui^} for which one of 
the following conditions holds: (a) there exists a 1 < j < p' such that x = Ti [j], for any x £ Ui 
and for any Ui £U; (b) x < Ti[l], for any x £ Ui and for any Ui £U; (c) Ti[p'] < x, for any 
X £ Ui and for any Ui £U; (d) there exists a 1 < / < p' — 1 such that Ti[j'] < x < Ti[j' + 1], 
for any x £ Ui and for any Ui £ U. For any set of related units U = {Ui^, Ui^, . . . , Ui^} we 
have that Usmc[^i/] = iy+i and Upr-e[^j/+i] = iy, for any I < y < z — 1. 

Concerning H" and H = HiHi . . . Hp'HpiHpi^iH' . Before this step all the elements in H were the 
original ones gathered in § 15.1.21 After the fourth step, the following will hold. 

{iv) The elements in H' and H" plus the elements in ffj[]HI[i] + 1 . . .u] and in ifj/[]H[[i'] + 1 . . . u], 
for any 1 < i < p' + 1 and 1 < i' < p', form the original set of elements that were in H before 
the fourth step. 

(v) We have that: (a) x < Ti[l], for any x £ -ffi[l . . . EI[1]]; (6) x > Ti\p'], for any x £ 
Hp'+i[l . . . M[p' + 1]]; (c) Ti[i-l] <x < Ti[i], for any x £ Hi[l . . . M[i\] and any 2 < i < p'; 
(d) X = Ti[i], for any x £ Hi[l . . . EI[i]] and any 1 < i < p' . 

(vi) (a) Let j = L[l] (j = L[p']), Uj is the rightmost unit such that x < Ti[l] (x > 

for any x £ Uj. (b) For any 2 < i < p', let j = L[i], Uj is the rightmost unit such that 
Ti[i — 1] < X < Ti[i] for any x £ Uj. (c) For any 2 < i < p' , let j = ]L[i], Uj is the rightmost 
unit such that x = Ti [i] for any x £ Uj. 

(vii) (a) N[l] (N[p']) is the number of x G 5[ such that x < Ti[l] (x > Ti\p']). (6) For any 2 < z < p', 
N[i] is the number of x G S'l such that Ti[i — 1] < x < Ti[i]. (c) For any 2 < i < p', N[i] is 
the number of x G 5[ such that x = Ti [i] . 

Let h = M[l], let i = 1 and let j = 1. We start the fourth step by scanning Bi. If Bi[i] < Ti[l], we 
increase h, and N[l] by 1, exchange B[i] with Hi[h] and increase i by 1. This sub-process goes 
on until one of the following two events happens: (a) h and are equal to u + l; (6) Bi[i] > Ti[l]. 
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If event (a) happens, we exchange the u elements currently in Hi with S'i[{j — l)u + l . . . ju]. Then, 
we set Upreb'] to L[l], I[Js„c[L[l]] to j and L[l] to j. After that, we set h and to and we 
increment j by 1. Finally, we go back to the scanning of Bi. Otherwise, if event (6) happens, we set 
h to ]HI[1] and we continue the scanning of Bi but with the following sub-process: if Bi[i] = Ti[l], 
we increase h, BI[1] and N[l] by 1, exchange B[i] with Hi[h] and increase i by 1. In its turn, this 
sub-process goes on until one of the following two events happens: (a') h and IHI[1] are equal to 
u + 1; {b') Bi[i\ > Ti[l]. Similarly to what we did for event (a), if event (a') happens, we exchange 
the u elements currently in Hi with S'i[{j — l)u + 1 . . ■ ju\. Then, we set Upreb] to L[l], UsuclLll]] 
to j and £.[1] to j. After that, we set h and ]HI[1] to and we increment j by 1. Finally, we go 
back to the scanning of Bi. Otherwise, if event (b') happens, we set h to EI[2] and we continue the 
scanning of Bi but with the following sub-process: if Bi[i] < Ti[2], we increase h, M[2] and N[2] by 
1, exchange B[i] with H2[h] and increase i by 1. We continue in this fashion, possibly passing to 
IHI[2], M[3], etc, until Bi is exhausted. Then, the whole process is applied to B2 from the beginning. 
When B2 is exhausted we pass to B^, B4 and so forth until each block is exhausted. 

Fifth. We start by exchanging H" and HiHi . . . Hp/Hp'Hp'^i. Let h = 1 and h' = u — M[l]. 
We exchange the elements in //i[IH[l] -|- 1 . . .u] with the ones in T2[h . . . h']. After that, we set 
h = h' , increment h' hy u — EI[1] and exchange the elements in lfi[IH[l] -|- 1 . . . u] with the ones 
in T2[h . . . h']. We proceed in this fashion until i^p'+i is done. Then we exchange the elements in 
H"H' with the rightmost \H"H'\ ones in T2. After that we execute the following process to "link" 
the HiS and HiS to their respective sets of related units. We start by setting Vpre\f' + 1] to L[l] 
and Us„c[L[l]] to /' + 1. Then we set VpreW + 2] to L[l] and U^^cfLfl]] to /' + 2. We proceed in 
this fashion until we set Vpre[f' + 2^ + 1] to l.\p' + 1] and Us„c[L[p' + I]] to f + 2p' + I. After 
that we execute the following process to bring each set of related units into a contiguous zone. Let 
i = [N[l]/nJ -|- 1, we exchange Hi with Ui; we swap the values in Upre[^] and \ipre[f' + l]i then the 
values in V^uM and Vsudf + 1]; we set Vpre[Vsuc[f' + 1]] = /' + 1 and Vsuc[Vpre[f' + 1]] = /' + !; 
finally, we set Us„c[Upre[^]] = i- Then, let j = Uprefi], we decrement i by 1, we exchange Uj 
with Ui] we swap the values in Upre[i] and Up^eb]) then the values b]; we set 

Upr-e[Usucb]] = 3 and Usnc[Upreb]] = finally, we set Upre[UsucH] = i and Usnc[Upre[i]] = i. We 
proceed in this fashion until the entire set of related units of Hi resides in ^^'[l . . . N[l] — 1HI[1]] [Hi 
now resides in 5^ [N[l] — EI[1] -|- 1 . . . N[l] — -|- u]). After that, we apply the same process to Hi, 
H2, H2 and so forth until every set of related units has been compacted into a contiguous zone. We 
end up with the sequence UiHiUiHi . . .U piHpiJJ piHpiJJ pU^iHpi^iH" H' , where each Ui contains 

the set of related units of Hi and each Ui the set of related units of Hi. Finally, we proceed to 
separate the elements of S[ from the ones residing in ffj[]HI[i] + \ . . .u] and i/j'[IHI[i'] -|- 1 . . . u], for 
any 1 < i < p' + \ and any 1 < i' < p'. Since the "intruders" were previously residing in T2, by 
the first and sixth steps in § 15.1.21 we know that any of them is less than any of the elements of 
S'l^. Therefore we can separate them with the stable partitioning algorithm in [6] and end up with 
the sequence RiRi . . . RpiRpiRpu^iH'" H" H' . Finally, we exchange the elements in the sequence 
H = H"'H"H' with the ones in T2, getting back in H its original (distinct) elements. 

Sixth. After the fifth step we are left with the sequence RiRi . . . Rp'Rp'Rp'^iH for which the 
following holds: (i) N[i] = \Ri\ and N[i'] = \Ri'\, for any 1 < i < p' + 1 and any 1 < i' < p'; (ii) 
for any x £ Ri, x < Ti[l]; (Hi) for any x G Rp'+i, x < Ti[p']; (iv) for any x £ Ri, x = Ti[i] 
{1 < i < p')] (v) for any x G Ri, with 1 < i < p' , Ti[i] < x < Ti[i + 1]. We begin by moving H 
before Ri with a sequence exchange. Since we do not need the p' pivots anymore, we put them 
back in their original positions in D[l . . .p] executing once again the process in the first step. Then, 
we allocate in the encoded memory an array M with p' + 1 entries of O {u) bits. Let i = 1 and 
M[l] = 1: ioi j = 2, . . . ,p' + 1 we increment i by Nb - 1] + ^[j - 1] and set Mb] = ^- After 
that, we execute a series oi p' + 1 recursive invocations of the procedure here in § 15.2.31 (not the 
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whole sorting algorithm). We start by sorting Ri = S'i\R[l] ... M [2] — 1] recursively with the same 
procedure in this section: we use Sj[M[l] . . . M[2] — 1] in place of 5[, D[l . . .p/t — 1] in place of D 
and so forth. After Ri is sorted, we swap H and Ri with a sequence exchange and proceed to sort 
i?2: we use 5;'[M[2] . . . M[3] - 1] in place of S'^, D[p/t + 1 . . . 2p/t - 1] in place of D and so forth. 
We proceed in this fashion until Rpi+i is sorted and H is located after it again. We do not need 
anything particularly complex to handle the recursion with O (1) words since there can be only 
0(1) nested invocations. 

Lemma 5. The aggregation phase requires 0{n) time, uses 0(1) auxiliary words and is stable. 
5.3 Final Phase 

The final phase has two main steps described in § 15.3.11 and § I5.3.2[ 
5.3.1 Sorting the main sequence S' 

After the aggregating phase we are left with S' = SiS'2 ■ ■ ■ S'^^ where, for any 1 < i < k, S[ is sorted 
and \S[\ = dlog^n. Let / = log^ n. We have three steps. 

First. We allocate in the encoded memory an array S'^ with | S'^ \ entries of one bit initially set to 
1. Then, we scan S'^. During the scan, as soon as we encounter a subsequence Si[i . . .i + f — 1] of 
equal elements (that is a subsequence of / consecutive equal elements) we set [i] , [i + 1] . . . S'^ [i + 
/ — 2] and §'^[i + / — 1] to 0. After that we use the partitioning algorithm in [6] to separate the 
elements of S'l with the corresponding entries in S' set to 1 from the ones with their entries set to 
(during the execution of the partitioning algorithm in [6], each time two elements are exchanged, the 
values of their entries in S' are exchanged too). After the partitioning we have that S[ = S'^Oi and 
the following conditions hold: (i) S'^ and Oi are still sorted {S'l was sorted and the partitioning is 
stable) ; {ii) the length of any maximal subsequence of consecutive equal elements of S'l is a multiple 
of /; {Hi) \Oi\ < df . Then, we merge Oi and 5*2 using the merging algorithm in [TU], obtaining a 
sorted sequence S2 with |52"| < 2d\og^ n. After that, we apply to S2 the same process we applied 
to S'l ending up with S'2' = S'202 where conditions (i), (ii) and (Hi) hold for S'2 and O2 too. We 
proceed in this fashion until each S"^ is done. 

Second. We have that S' = S"0. The following conditions hold: (i) O is sorted and \0\ < df; 
(ii) the length of any maximal subsequence of consecutive equal elements of S'{ is a multiple of 
/ (and so | S" \ is a multiple of / too) . Let s = \S" \ / f and let us divide S" into s subsequences 
F1F2 . . . Fs-iFs with \Fi\ = f. We allocate in the encoded memory two arrays Spj,g and S'/^ci each 
one with s entries of Q (logn) bits. We also allocate an array C with r (= n'') entries of Q (logn) 
bits, each one initialized to 0. Then, for each Fi from the rightmost to the leftmost one, we do the 
following: Let v = Fi[l]. If C[v] = 0, we set Sg^cH = and C[v] = i. Otherwise, if C[v] ^ 0, we set 
S'4cH = CM, %re[C[v]] = i and CM = i. 

Third. We scan C and find the leftmost entry not equal to 0, let it be i. Let j = C[z], we 
exchange Fi with Fj and do the following: (i) we swap the values in Spre[l] and §pVe[j]) then the 
values in S'^Jl] and (ii) we set ^'preWLAjW = J and S"„JSpVeb1] = J; (^^^) finally, we set 

Spre&c[l]] = 1 and S';„c[Spre[l]] = 1- Then, let j' = S';„Jl]. We exchange F2 and Fj, and then 
we make similar adjustments to their entries in S^'^g and We proceed in this fashion until we 
exhaust the linked list associated with the ith entry of C. After that we continue to scan C, find 
the leftmost non-zero entry i' > i and process its associated list in the same way. At the end of the 
process, the Fi^s have been permuted in sorted stable order. Now both S" and O are sorted. We 
merge them with the merging algorithm in [10] and S' is finally sorted. 
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5.3.2 Taking care of the encoded memory and the zones 

In the last main step of the final phase we sort all those zones that have been built in the preliminary 
phase. We scan G, let u G G be the ith. element accessed, we exchange Ti[i], T2[i] and T-^li] with 
Ji[u], J2[u] and Jslu], respectively. We swap S" and H {H has been moved after S" at the end 
of the aggregating phase). W, G, D, P, Ti, T2, T3 and H have O (d) = O {rf) elements and we 
can sort them using the mergesort in |10] . The obtained sequence can now be merged with S" 
using the merging algorithm in [10]. Ji, J2, J3 and V have O(n^) elements. We sort them using 
the mergesort in [10]. Finally, M' and M" have 0(n/logn) elements. We sort them with the 
mergesort in [TO] and we are done. 

Lemma 6. The final phase requires 0{n) time, uses 0(1) auxiliary words and is stable. 
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